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The magnetization curve of the two-dimensional spin-1/2 J1-J2 Heisenberg model is investigated 
by using the Chern-Simons theory under a uniform mean-field approximation. We find that the 
magnetization curve is monotonically increasing for J2/J1 < 0.267949, where the system under zero 
external field is in the antiferromagnetic Neel phase. For larger ratios of J2/J1, various plateaus will 
appear in the magnetization curve. In particular, in the disordered phase, our result supports the 
existence of the M/Msat = 1/2 plateau and predicts a new plateau at M/Msat = 1/3. By identifying 
the onset ratio J2/J1 for the appearance of the 1/2-plateau with the boundary between the Neel and 
the spin-disordered phases in zero field, we can determine this phase boundary accurately by this 
mean-field calculation. Verification of these interesting results would indicate a strong connection 
between the frustrated antiferromagnetic system and the quantum Hall system. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.50.Ee, 73.20.Dx 



Due to quantum and frustration effects, rich physics 
can appear in thej-|frustrated quantum spin systems at 
zero external fieldJJ Exciting behavior was also observed 
recently in several cases with an external magnetic field. 
For instance, the recently discovered two-dimensional 
(2D) S = 1/2 spin-gap material SrCu2 (603)2, ^hich 
can be described by the Shastry-Sutherland modelp ex- 
hibits several plateaus at M/Mgat = 1/3, 1/4, and 1/8 
in its magnetization curve, where M (Msat — 1/2) is 
the (saturating) magnetization per site.H The origin of 
these plateaus and the nature o£±hfi_corresponding spin 
states are under intense debatcaD'EllII Recently, by map- 
ping onto spinless fermions carrying one quantum of sta- 
tistical flux and under a mean-field approximation, Mis- 
guich et al. show that the original spin model can be 
related to a generalized Hofstadter problem, where the 
spin excitation gaps that produce the observed magneti- 
zation plateaus arise from some of the Landau level gaps 
in the integer quantum Hall effect for the fermions on a 
lattice.El For realistic values of the exchange constants, 
they obtain an excellent quantitative fit to the observed 
magnetization curve, which demonstrates the success of 
their approach. 

Another prototype of a realistic frustrated two- 
dimensional system, which has been recently real- 
ized experimentally in Li2VOSi04 and Li2VOGe04 
compounds,cl is the so-called J1-J2 Heisenberg model 
with the Hamiltonian 



- ^ Jij Si ■ Sj - S-, 



(1) 



where the exchange couplings Jy are equal to Ji when 
i,j are nearest neighbors on the square lattice; Jij 
are equal to J2 when i,j are connected by a diago- 
nal bond. The external magnetic field B is applied 
along the z-axis. Both couplings are antiferromagnetic. 
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i.e. Ji,2 > 0, and the spins Si = 1/2. This model 
has 
years. I 

J2 / Ji where the frustration is weak, the system ex- 
hibits Neel ordering described by a wave vector (tTjTt). 
When J2/J1 is large enough, the ground state is domi- 
nated by the next-nearest-neighbor interactions and has 
a colHnear order described by (7r,0) or (0, tt). There is 
also a general consensus on the disappearance of the mag- 
netic ordering at 0.38 < J2/J1 < 0.6, while the identi- 
fication of thi; gramid^ ^atje is still a subject of much 
controversy.HByyy.By 

Just like the case of the Shastry-Sutherland model, 

curve are predicted 
while the situation is 
even more controversial. Strong evidence for a plateau 
at M/Msat = 1/2 in the region 0.5 < -h/^r^ 0.65 
has been recently reported by Hanecker et a/.c3E2l How- 
ever, Fledderjohann and Miittertd do not find a plateau 
at M/Msat = 1/2 in the region of the quantum disor- 
der phase; but instead they find some indications for a 
plateau-like structure at M/Mgat = 2/3. Most of the 
previous studies on the plateau are obtained by the nu- 
merical calculations on a small clusters (the typical num- 
ber of lattice sites in these works is about 6x6 = 36). 
As discussed in Ref. |2], the plateau structures in the 
magnetization curve can depend sensitively on the sys- 
tem size. Therefore, these results may be plagued with 
the finite-size effects. For example, some of the predicted 
plateaus may be an artifact of the special lattice geome- 
try, and the boundary conditions used may frustrate the 
order which tends to develop. Thus the precise deter- 
mination of the positions and widths of the plateaus in 
the frustrated Heisenberg model is indeed a very delicate 
problem, and a better theoretical understanding of the 
magnetic order and of the mechanisms which create the 
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plateaus is needed. 

To avoid the possible finite-size effects, we apply 
the Cheijn-Simons (CS) theory for the magnetization 
plateausia to the 2D spin-1/2 J1-J2 Heisenberg model. 
Because of its success for the Shastry-Sutherland model, 
this approach should give us reasonable results in the 
present case. We find that the magnetization curve is 
monotonically increasing for J2/J1 < 0.267949, where 
the system at zero field is in the antiferromagnetic Neel 
phase. Besides, various plateaus will appear in the the 
magnetization curve for larger ratios of Jij J\. In partic- 
ular, in the disordered phase, our result supports the ex- 
istence of the M /Msat = 1/2 plateau and predicts a new 
plateau at M/Mgat = 1/3. Furthermore, we note that, by 
identifying the onset value of J2/ Ji for the appearance of 
the 1/2-plateau with the critical value of the phase tran- 
sition between the Neel and the spin-disordered phases 
in zero field, this phase boundary, \5fhjch was determined 
earlier by heavy numerical meansjljll3 can be reproduced 
accurately by this mean-field calculation. 

The Hamiltonian in Eq. (|^) can be rewritten as 

H = Hxy + Hz — B ^ , 

i 

H.y = \Y. {S^Sr + S+Sr) , (2) 
Hz = Jij S^Sj, 

where H^y and H^ are the spin-flip and the Ising parts 
of the Hamiltonian at zero field. According to Ref. ||, 
H,^y and H^ are treated in different ways. First, the 
Ising part is approximated by a simple uniform mean- 
field decoupling, 

Hz ~ 4( Ji -I- .h)M St - 2( Ji + J2)M^N, (3) 

i 

where N is the number of lattice sites. Thus Hz gives a 
contribution to the total energy with a simple dependence 
on the magnetization. Second, for H^y, one can exactly 
map the spin operators to spinless fermions attached with 
a fiux tube carrpfirig .one flux quantum of statistical CS 
magnetic fieldE^'Oa where 5'f -I- 1/2 corresponds to the 
occupation number rii of site i. Under a mean-field treat- 
ment such that the flux tubes are smeared out into a 
uniform background magnetic field, the fiux per square 
plaquette 4> is then tied to the density of fermions and 
thus to the magnetization M of the spin system : 

l^{n)^M^\. (4) 

Because of this fiux, each energy band splits to subbands 
with a complicated structure. Therefore, the present spp. 
system can be identified with a Hofstadter problemCJ 
for fermions moving on a square lattice witlij.nearest- 
neighbor and next-nearest-neighbor hoppings.EZrE3 For a 



given M [or (/) = 2-k{M + 1/2)], the mean-field ground 
state is obtained by filling the lowest energy subbands 
with fermions until their density satisfies (n) — (/)/27r. 
The one-body problem from H^y can be straightfor- 
wardly analyzed for rational values of 0/27r. For 0/27r 
p/q {p and q are mutually prime integers), there are q 
subbands and the ground state is the Slater determinant 
with the lowest p subbands being completely filled. The 
energy of the filled subbands leads to another contribu- 
tion to the total energy. Hence the total energy per site 
E{M) as a function of the magnetization becomes 

= ^ E E4"^+2(^i + ^2)M^ (5) 

Q = l,....p £ 

Here el"-* is the eigenenergy of the a-th subband with the 

wavevector k being restricted to the magnetic Brillouin 
zone. The magnetization can be obtained as a function 
of B by minimizing E{M) — BM. It is clear from Eq. (|^) 
that, without the contribution from H^y, the magnetiza- 
tion M will be linearly proportional to B, and there is 
no magnetization plateau. Therefore, the appearance of 
magnetization plateaus is related to certain features of 
the Hofstadter spectrum. 

The Hofstadter diagrams for J2/Ji~0.2 and 0.3 are 
shown in Fig. ^ where the lower bold line marks the 
Fermi level (highest occupied state) and the upper one 
marks the lowest unoccupied level. A jump of the Fermi 
energy as a function of M in Fig. 1(b) leads to discon- 
tinuity of the slope of the function E{M). These jumps 
for various M give rise to plateaus in the magnetiza- 
tion curve. They are closely related to the occurrence 
of band-crossing when the value of J2I Ji is varied. For 
example, in Fig. 1(a) before the upper two subbands at 
<l)/2iT = 2/3 touch at J2/J1 = 0.267949, the Fermi en- 
ergy is continuous and there is no magnetization plateau 
at Af/Afsat = 1/3 (see Fig. |). However, in Fig. 1(b), 
the "pockets" enclosed by the bold lines are separated 
at 0/27r = 2/3. It can be seen that the Fermi level to 
the right of the contact point is below (above) the band 
gap before (after) band-crossing. As discussed earlier, 
the Fermi level marks the position of the p-th subband 
[see Eq. (|^)]. Therefore, apparently a subband associ- 
ated with some fiux slightly larger than 2/3 is shifted 
above the energy gap after band-crossing. This shift of 
a fine subband due to the crossing of the broader sub- 
bands at (/)/27r = 2/3 was studied earlier in the con- 
text of the Hofstadter spectrum.ca It is closely related 
to the jump of the iii|tepe|r-valued Hall conductances of 
the broader subbands .E3c3 The emergence of the magne- 
tization plateaus for the spin system thus has an inter- 
esting connection with the change of the integer-valued 
Hall conductances induced by band-crossing. 

To justify the present approach, it is important to 
check whether the plateau states are robust against fiuc- 
tuations around the mean-field solutions. It was showed 
that the Gaussian fluctuations of the CS, gauge field are 
massless only when the TKNN integeiO describing the 
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FIG. 1: Hofstadter spectra for a square lattice with nearest- 
neighbor and next-nearest- neighbor hoppings for J2/J1 = 0.2 
(a) and 0.3 (b). Vertical lines mark the energy bands as 
a function of the statistical flux 0/27r per square plaque- 
tte. Total Hall conductances ctt above the Fermi level from 
M/Msat = {(l>/2n = 1/2) to M/M^at = 1 (0/27r = 1) are in- 
dicated. The Hall conductance at 4)/2tt = 2/3 changes from 
1 to — 2 ffllpgn the upper two subbands touch at J2/J1 
0.267949 .Bt 



quantized Hall coefficient of-.the fermions on the frus- 
trated lattice becomes unity£3 In that case, the Gaussian 
fluctuations induce instabilities for the mean-fleld ground 
states. We have computed the TKNN integers numeri- 
cally (for example, see Fig. |^) and found the Gaussian 
fluctuations to be massive. Therefore, the plateaus are 
not destroyed by quantum fluctuations. 

Magnetization curves for various J2/J1 ratios are 
shown in Fig. |2[ We note that the saturation field Ssat 
can be computed exactly by identifying the energy Ef of 
the fully polarized state with the (exact) minimum en- 
ergy -E™" of the states with one spin flipped, Ej = E'™". 
Thus B,at/Ji = 4 for J2/J1 < 1/2; B,at/Ji = 2 + 4J2/J1 
for J2/J1 > 1/2. Our flndings agree with these exact 
results near the saturation. 

In the Neel phase, it is expected that the spins cant 




FIG. 2: Magnetization curves for the J1-J2 Heisenberg model 
calculated using uniform CS mean-field. The curves from left 
to right are for J2/J1 — 0,0.2,0.3,0.4,0.5, and 0.7, respec- 
tively. 



gradually from the antiparaUel configuration toward the 
parallel configuration until the magnetization saturates 
at the saturation field Bgat ■ The magnetization curve ob- 
tained from the present approach is consistent with this 
expectation: it is featureless all the way to full saturation 
when J2/J1 is small (see the curves for J2/J1 — and 
0.2 in Fig. Upon increasing J2/J1, plateaus emerge 
and the magnetization curves become more complex. In 
particular, because of the band-crossing at (j)/2TT = 2/3 
when J2/J1 = 0.267949 (see Table I of Ref. ||), a plateau 
at M/Msat = 1/3 is found [see Eq. (|)]. This is an unex- 
pected result, especially for 0.267949 < J2/J1 < 0.38 
where the ground state in the absence of an external 
field is |in_t|^_CvIe£^-phase. In the previous finite-size 
studies,BEaoEl'BB there is no indication for the ap- 
pearance of this plateau. However, the system sizes and 
the boundary conditions they used forbid the appearance 
of the M/Msnt = 1/3 plateau, therefore the possibility of 
this plateau cannot be ruled out. Furthermore, when 
J2/J1 — 0.382683, another band-crossing in the Hofs- 
tadter spectrum occurs at 0/27r = 3/4 (see Table I of 
Ref. |2^), and a plateau at M/Ms^t = 1/2 ensued. It 
is interesting to find that the critical value for the ap- 
pearance of the M/Msat = 1/2 plateau agrees quite well 
with the critical point of the quantum phase transition 
at zerOrfijdd between the Neel and the quantum disorder 
phases .Mr 'I This reinforces our confidence on the present 
CS mean-fleld approach. As mentioned before, while the 
appearance of a plateau in the quantum disorder phase 
had beeBiipredicted, the value of the plateau is still under 
debate.c2lLifE3 The controversy may come from the subtle 
flnite-size effects in their investigations. Since the present 
CS theory is free from the flnite-size effects, we give a 
strong support for the existence of the M/Msat = 1/2 
plateau. 
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More complex structures in the magnetization curves 
appear when J2/J1 is further increased. For example, 
when J2/J1 = 1/2, a series of plateaus at M/ A/gat = 
n/{n + 2) is found, which corresponds to the band- 
crossing at some magic numbers 0/27r — (n + l)/(n + 2). 
Irregular plateau structures are found for even higher 
values of J2/J1 (see the J2/J1 = 0.7 curve in Fig. ||). 
This behavior is quite similar to the case of the trian- 
gular lattice, where many plateaus are ppedicted under 
the uniform CS mean-field approximation.Ll In the case of 
the triangular lattice, it is shown that, when the nonuni- 
form mean- field solutions are used, only the main plateau 
{M/Msat = 1/3 in that case) survives and other mini- 
plateaus disappear. We believe that the same situation 
will happen in the present study of the J1-J2 model. That 
is, when going beyond the uniform mean-field approxi- 
mation, many of the mini-plateaus in the magnetization 
process may disappear and only the main plateaus with 
simple fractions of M /Msa± survive. Moreover, as shown 
in Fig. 1^, the spin gap (the plateau at M/Mgat = 0) opens 
at J2/J1 = 0, instead of opening at the correct critical 
coupling J2/J1 ~ 0.38. This feature comes from the fact 
that the Neel state of the square-lattice antiferromagnet 
is not correctly described by the present uniform mean- 
field approximation. 

In conclusion, the magnetization curve of the 2D spin- 
1/2 J1-J2 Heisenberg model is studied by using the 



Chern-Simons theory under a uniform mean-field approx- 
imation. In the disordered phase, our result supports 
the existence of the M/Mgat = 1/2 plateau and pre- 
dicts a new plateau at M/Mgut — 1/3. Moreover, various 
plateaus appear in the magnetization curve both in the 
disordered and the collinear phases, which could be an 
artifact of our mean field approximation. More work is 
needed to be conclusive at this point. We note that it is 
experimentally accessible to confirm our results. As men- 
tioned before, the 2D spin-1/2 J1-J2 Heisenberg model 
with J2/J1 — 1 has been recently realized experimen- 
tally in Li2VOSi04 ( Ji + -^2 rf 8-2 K) and Li2VOGe04 
( Ji + J2 — 8.2 K) compounds.B If the g-factor is taken to 
be 2, the corresponding saturation fields will be approx- 
imately 18 T for Li2VOSi04 and 13 T for Li2VOGe04. 
In both cases, these values of the magnetic fields can 
be reached experimentally. Thus the full magnetization 
curve could be mapped out to test our results. Veri- 
fication of the magnetization plateaus would indicate a 
strong connection between the frustrated antiferromag- 
netic system and the quantum Hall system. 
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